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Shear-induced anisotropic decay of correlations in hard-sphere col- 
loidal glasses 
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Abstract - Spatial correlations of microscopic fluctuations are investigated via real-space experi- 
ments and computer simulations of colloidal glasses under steady shear. It is shown that while the 
distribution of one-particle fluctuations is always isotropic regardless of the relative importance 
of shear as compared to thermal fluctuations, their spatial correlations show a marked sensitivity 
to the competition between shear-induced and thermally activated relaxation. Correlations are 
isotropic in the thermally dominated regime, but develop strong anisotropy as shear dominates 
the dynamics of microscopic fluctuations. We discuss the relevance of this observation for a better 
understanding of flow heterogeneity in sheared amorphous solids. 
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Introduction. — The nature of microscopic fluctuat- 
ions in the relaxation and flow of glasses is a topic of 
great current interest because of their importance to a 
wide range of materials including foams, emulsions, granu- 
lates, and colloidal and molecular glasses. Thermal glasses 
such as molecular and colloidal glasses exhibit structural 
relaxation due to thermally induced rearrangements [TJ|3] , 
while athermal amorphous materials such as granulates 
and foams relax typically by external forces such as ap- 
plied shear that drive their flow Hip]. It has been sug- 
gested that the applied shear acts as an effective temper- 
ature that drives microscopic fluctuations in an isotropic 
way [6-8 . Such a concept is attractive since it opens a 
possibility of the generalized fluctuation dissipation rela- 
tion under a nonequilibrium situation. This idea is also 
motivated by the jamming phase diagram p\ that depicts 
applied shear stress as an alternative way to fluidize the 
glass, in an attempt to unify amorphous materials in a 
common framework. Supported by simulations of sheared 
glasses [6 p], t he current models of plasticity of amorphous 
materials |1Q[[TT] assume that the rate of activated events 
is controlled by an isotropic effective temperature that is 
dependent on the shear rate. While these models have 
been quite successful in describing many aspects regard- 
ing the plastic response of amorphous solids to external 
load such as simple shear, many other issues such as flow 



heterogeneity and shear banding still remain under de- 
. Some additional ingredient seems, therefore, 
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to be necessary for a full description of the deformation 
behavior of amorphous solids. 

Here, we focus on the issue of anisotropy. Since the 
applied shear imposes a directionality of motion, this 
should reflect itself in the nature of microscopic fluctu- 
ations and/or their correlations. In fact, some evidence 
of direction-dependent correlations comes from quasistatic 
computer simulations of two-dimensional Lennard- Jones 
glasses at zero 15|[l6 and finite [17] temperatures. Direct 
observation of these fluctuations in molecular glasses, how- 
ever, is prohibited by the small molecular length scales. 
We use real-space experiments and molecular dynamics 
simulations of colloidal hard spheres to investigate the ef- 
fect of shear on the microscopic fluctuations. Hard-sphere 
colloidal glasses provide good model systems for molec- 
ular glasses, with the advantage that the motion of the 
individual particles can be followed directly in real-space 
and time using confocal microscopy. The direct imag- 
ing of displacement fluctuations allows us to identify a 
novel kind of correlation with stress-dependent anisotropic 
scaling. We determine correlations by following fluctu- 
ations of non-affine displacements directly in real space, 
and we investigate their direction-dependence. By vary- 
ing the applied shear rate, we probe two different relax- 
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at ion regimes: the slowly driven "viscous" regime where 
the shear rate 7 is smaller than the inverse structural re- 
laxation time, r, so that relaxation occurs predominantly 
by thermal fluctuations, and the shear-dominated regime 
where 7 > r _1 , so that relaxation is driven by the applied 
shear. In agreement with the idea of a scalar effective 
temperature, we find that the single particle displacement 
fluctuations are always isotropic, even when shear domi- 
nates the relaxation. At the same time, however, we ob- 
serve a striking difference in the spatial correlations of fluc- 
tuations: while in the slowly driven regime, correlations 
exhibit isotropic decay, i.e. with a direction-independent 
scaling exponent, in the shear-driven regime, this decay 
becomes anisotropic, being slowest in the flow direction. 
We argue that such stress-dependent anisotropic correla- 
tions can cause the ubiquitous shear banding instability 
of amorphous materials. These results reveal a new kind 
of criticality exhibited by amorphous solids on the verge 
of flow. The particular feature of this stress-induced crit- 
icality is the stress-dependent anisotropy of the scaling of 
correlations. 

Experiments. — Our colloidal glass consists of ster- 
ically stabilized polymethylmethacrylate (PMMA) parti- 
cles suspended in a mixture of Cycloheptyl Bromide and 
Cis-Decalin that matches both the density and refractive 
index of the particles. The particles have a diameter of 
<j = 1.3/im, with a polydisperity of 7 %. The particle vol- 
ume fraction is <j) ~ 0.6, well within the glassy state (note 
that </> g iass ~ 0.58) 18 . The suspension is loaded in a cell 
between two parallel plates at a distance of L z = 65/im. 
Shear is applied by moving the top plate in the positive 
x— direction using a piezoelectric translation stage. Con- 
focal microscopy is used to image the motion of the indi- 
vidual particles in a 108/im x 108/im x 65/im volume, and 
to determine their positions in three dimensions with an 
accuracy of 0.03/im in the horizontal, and 0.05/im in the 
vertical direction [2]. The motion of ~ 2 x 10 5 particles is 
then tracked during a time interval of 25 min by acquiring 
image stacks every 60s. All measurements are taken in the 
steady state regime, after the glass has been sheared to a 
strain of 7 ~ 1. 

As a basic quantity, we determine non-affine displace- 
ments as follows 
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For each particle, we follow nearest 
neighbors in time and determine the best affine tensor e 
that transforms the nearest neighbor vectors, di = ri — Vq, 
over the time interval St. This is done by minimizing 
D 2 = (l/n)Y^ =1 (di(t + St) - (I + e) ■ di(*)) 2 , where I is 
the identity matrix. D 2 reflects the mean-square deviation 
from a local affine deformation, and is an excellent mea- 
sure of local plasticity [l0]. The non-affine displacements 
of the particles are determined both (i) via subtraction of 
the local flow, dr na = r(t + St) - r(t) - f dt'u(t', r{t')), 
where u(t,r(t)) is a coarse-grained displacement field [l9] 
and (ii) as dr na = r{t + St) — r(t) — e • r(t). It is assumed 
in (ii) that the coordinate center is at rest. We have veri- 
fied that these two definitions give identical results as long 
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Fig. 1: Fluctuations of single particle displacements around 
the mean flow (Ax, blue symbols) and in the shear- gradient 
direction (Az, red symbols). Circles and squares indicate 
the thermal (Pe* = rj = 0.8), and shear-dominated regime 
(Pe* = 4.4), respectively. All distributions overlap, demon- 
strating perfect isotropy. 



as the shear is homogeneous across the channel [20]. In 
the general case of heterogeneous shear, however, the first 
definition is used. 

To separate the effects of shear and temperature, we dis- 
tinguish applied shear rates by the modified Peclet num- 
ber [21] Pe* = 7T, the product of the applied shear rate 
7 and the structural relaxation time r of the glass. When 
Pe* < 1, the applied shear rate 7 < r _1 , and the re- 
laxation is dominated by thermal fluctuations. However, 
when Pe* > 1, 7 > r _1 , and relaxation is dominated by 
the applied shear. We define the structural relaxation time 
of the quiescent glass as the time that the mean square 
displacement of particles exceeds the particle radius. For 
our system, r ~ 2 x 10 4 s, a factor of 5 x 10 4 larger than 
the Brownian diffusion time t# = 0.4s in the dilute limit. 
We apply shear rates between 7 = 1.5 x 10 _5 s _1 and 
2.2 x 10 _4 s _1 that correspond to Peclet numbers between 
0.3 and 4.4, smaller and larger than one. This allows us to 
investigate the effect of temperature versus shear on the 
correlations of microscopic fluctuations. 

Experimental results and discussion. — In agree- 
ment with the idea of an effective temperature, the dis- 
placement fluctuations of the individual particles are al- 
ways isotropic, as demonstrated in Fig. [TJ The figure 
shows the fluctuations of non-affine displacements in the 
flow (x) and the shear-gradient direction (z) for both 
Pe* < 1 and Pe* > 1. Clearly, the distribution of dis- 
placements in both directions are identical, and symmet- 
ric with regard to forward and backward jumps. These 
properties hold regardless of the Peclet number, i.e. both 
where thermal fluctuations determine the dynamics, and 
where shear is the dominant factor. We conclude that 
there is no shear-induced anisotropy in the fluctuations 
of the individual displacements, indeed favoring an "effec- 
tive temperature" picture, in which the applied shear is 
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Fig. 2: Pair correlation function g(r) at Pe* > 1, where shear 
dominates the relaxation, (a) Two-dimensional map resolving 
g(r) in the plane perpendicular to the shear axis. Dashed circles 
demonstrate isotropic first and second nearest neighbor shells, 
(b) Angle-resolved pair correlation function as a function of dis- 
tance along the 45° (axis of extension) and 135° direction (axis 
of compression). Bins of width 7r/18 around these directions 
were used for averaging. 



accounted for by an isotropic effective noise. 

Regarding the structure of the glass, some anisotropy 
is expected due to the presence of a finite shear stress 



22 . This anisotropy is, however, quite small in the range 
of volume fractions and shear rates of interest to this 
work. To demonstrate this, we determine the direction- 
dependent pair correlation function g(r) that indicates 
the probability of finding two particles separated by r. 
A two-dimensional density map that resolves g(r) in the 
plane perpendicular to the shear axis is shown in Fig. [2^i. 
The pair correlation function appears closely circular sym- 
metric, indicating that the structure remains -within 
our experimental resolution- essentially isotropic. The 
slight deviation from circular symmetry is caused by the 
lower instrument resolution in the vertical direction; we 
checked this by imaging the glass without applied shear, 
where we obtain a very similar g(r). Nevertheless, some 
anisotropy is expected with respect to the two diagonal di- 
rections [22]; this anisotropy is, however, very small. We 
resolve g{r) along 45° and 135° to the horizontal, and plot 
the angle- specific g(r) in Fig.[2]3. The close overlap of the 
two curves shows that the anisotropy is of the order of our 
experimental accuracy, and can barely be resolved. We 
note that this is in contrast to measurements performed 
at much higher Pe where the applied shear can have a 
pronounced effect on the structure, and the pair correla- 
tion becomes clearly anisotropic [23]. At the moderate 
Pe investigated here, however, both the structure and the 
displacement fluctuations remain essentially isotropic and 
barely any signature of the applied shear is observable. 

In contrast to the static correlation function and single- 
particle displacements, a surprisingly strong signature of 
the applied shear is observed in the two-point correlations 
of displacements. To demonstrate this, we determine the 




Fig. 3: Anisotropy of spatial correlations. Correlation func- 
tion C D 2 resolved along the horizontal (red circles) and the 
vertical directions (black squares and blue diamonds), in the 
thermal (a) and shear dominated regimes (b). Bins of width 
7r/18 around the horizontal and vertical directions were used 
for averaging. C D i shows isotropic algebraic decay in the ther- 
mal (a), and anisotropic decay in the shear dominated regime 
(b). Insets show same data on linear scale. 



spatial correlation function [24] 



C c2 (Ar) 



(D 2 (r + Ar)D 2 (r)) - (D 2 (r)) 2 
P 2 (r) 2 ) - <£ 2 (r)> 2 



(1) 



which correlates non-affine displacements at locations sep- 
arated by the vector Ar = (5x,5y,5z). C^2 correlates 
fluctuations of the plastic activity, and is analogous to cor- 
relation functions defined in critical phenomena for equi- 
librium fluctuations. We have recently observed that Cx>2 
exhibits a power-law decay over the entire system size, 
demonstrating that plastic flow is organized in a hierarchi- 
cal way [24]. Here, we elucidate the direction-dependence 
of correlations by resolving C^2 along the horizontal and 
vertical directions. To do so, we average Cj^i in angular 
wedges of width tt/18 and plot the angle-specific corre- 
lation function in Fig. [3] A striking difference between 
the thermal and the strongly driven regime is observed. 
While in the thermal regime, correlations decay isotropi- 
cally with the same power law Cd 2 ~ r~ a with exponent 



p-3 



V. Chikkadi et al. 



a ~ 1.3 (Fig. [3^l), in the stress-dominated regime, the de- 
cay becomes anisotropic: correlations decay faster in the 
vertical, and slower in the horizontal direction (Fig. [3)3). 
These results suggest that the application of shear on col- 
loidal glasses leads to the manifestation of new correlations 
15 with shear-rate dependent anisotropic decay. 



Simulations. — In our experiments, flow is imposed 
at the boundaries, and the glass flows inhomogeneously, 
i.e. develops shear bands in the shear-dominated regime; 
the question is then whether these anisotropic correlations 
are the result of the boundary conditions. To investigate 
this, we use computer simulations, where the boundary 
conditions can be varied more easily than in the experi- 
ment. We perform event driven MD simulations of a poly- 
disperse (11%) hard sphere system in three dimensions. 
The quiescent properties of this system have been studied 
extensively in 18,25]. The temperature is fixed at T = 1 
via velocity rescaling (in our simulations, particle mass 
and diameter and the Boltzmann constant are set to unity, 
thus setting the unit of time and energy also to 1). In the 
simulations, N = 30000 particles are placed in a random 
configuration between two walls, separated by 30 mean 
particle diameters. A constant shear is applied by mov- 
ing the walls with velocities ±Z7 wa ii in indirection. We 
perform simulations in the thermal and shear-dominated 
regimes by keeping the shear rate constant (7 = 10 -4 ) but 
varying the relaxation time, r, via a change of the volume 
fraction from <j> = 0.57 (r « 10 3 -» Pe* « 0.1) to a value 
in the glassy state, <j> = 0.60 (r > 10 5 Pe* > 10). In ad- 
dition to the wall-driven flow, we have also performed sim- 
ulations using Lees-Edwards boundary conditions, where 
the top of the simulation cell is identified with the bottom 
displaced with a constant horizontal velocity 26 . While 
the wall-driven flow better mimics the experimental situ- 
ation, the latter has the advantage of being free of wall 
effects, thereby revealing the bulk response of the glass 
27. 

The results of these simulations are shown in Fig. [4] 
Interestingly, but not unexpectedly, due to rather low 
stresses present in these studies, only a slight anisotropy is 
observed in the structure (Fig. |4^l). Furthermore, in very 
good agreement with our experiments as well as previous 
reports [29], the fluctuations of the individual displace- 
ments exhibit essentially isotropic behavior (Fig. |4Jd). We 
observe a gaussian central part and a perfect exponen- 
tial decay at large displacements; the universal character 
of these features has been the subject of recent studies 
(see, e.g., [28] and references therein). In contrast, strong 
anisotropy is observed in the correlations of displace- 
ments when the Peclet number increases. Correlations 
decay isotropically in the thermal regime (Fig. [4^), while 
in the shear-dominated regime, the decay is anisotropic 
(Fig. |4]i). This is true for both wall-driven as well as 
Lees-Edwards-boundary conditions, indicating that these 
anisotropic correlations are rather a bulk property of the 
driven glass. This anisotropy seems to be a robust feature 
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Fig. 4: Results of wall-driven computer simulations in the ther- 
mal and shear-dominated regimes. The shear rate is 7 = 10 -4 , 
corresponding to Peclet numbers Pe* = 7-7 « 0.1 for = 0.57 
and Pe* > 10 for <j) > 0.60. (a) Pair distribution function in 
the shear plane (spanned by the flow and the shear-gradient 
vectors) along lines of 45° (axis of extension), 90° (flow gradi- 
ent direction) and 135° (compression axis) with respect to the 
flow direction. Results for the flow direction (0°) are indistin- 
guishable from those of the flow gradient direction. Error bars 
give standard deviation over 13 independent simulation runs, 
(b) Distribution of single particle displacements along the flow 
(Ax, after subtraction of the flow contribution) and flow gra- 
dient direction (Az). The time interval for the measurement 
of Ax and Az corresponds to 1% strain. Dashed (solid) lines 
indicate gaussian (exponential) fits to small (large) displace- 
ments 28 . (c) and (d) Spatial correlations of plastic activity, 
Eq. |l]), along the directions parallel and perpendicular to the 
flow. 
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of shear-dominated flow: while the spatial decay of corre- 
lations deviate from a power law, possibly due to the differ- 
ent boundary conditions, the change of symmetry persists, 
indicating that the anisotropy is a fundamental property 
of shear-induced relaxation. 



Flow instabilities. — This anisotropic decay can 
have important consequences for the mechanical stability 
of amorphous materials. We hypothesize that it can lead 
to the ubiquitous shear banding instability: the longer 
range of correlations in the horizontal direction leads to 
enhanced coupling of non-affine displacements in this di- 
rection, which in turn can lead to highly correlated flow 
that spans the sample in the form of a shear band. At the 
same time, the faster decay of correlations in the vertical 
direction restricts the coupling of particle displacements, 
and constrains the shear band to a rather small thickness 
(see Figs, ^jp and |4]i). To test this idea, in the exper- 
iment, we investigated whether the anisotropy of corre- 
lations persist during the initial transient regime, where 
the flow is still uniform. Indeed, we find that already in 
this transient regime before any shear banding manifests, 
correlations are strongly anisotropic, suggesting that it is 
this anisotropic nature of correlations that causes the well- 
known shear banding instability. 



Conclusion. — Our experiments and computer sim- 
ulations reveal a surprising stress dependence of micro- 
scopic fluctuations in glasses. Correlations of non-afflne 
displacements are anisotropic, and exhibit unusual, stress- 
dependent decay. This decay is direction-dependent when 
shear dominates the relaxation, and becomes isotropic 
when thermal relaxation prevails. These results have im- 
portant implications for models and theories of the plastic 
deformation of amorphous materials. While unifying con- 
cepts such as effective temperature p] and the jamming 
phase diagram 9 have been quite successful in describing 
many aspects of the response of driven glassy systems, our 
results show that a more refined description should also 
account for the shear-induced anisotropy in spatial cor- 
relations. This shear-induced anisotropy suggests a new 
route for understanding the origin of flow heterogeneities 
in a wide range of amorphous materials. Because at her- 
nial granular materials are always dominated by the ap- 
plied shear, they show always shear localization [5 15 



while thermal amorphous materials show shear instabil- 
ities only in the glassy phase but flow uniformly in the 



supercooled state 12 13,24 
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